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Abstract 

We present a preliminary search for the flavor-changing neutral current decay B + — ► K + vV using 
56 million T(4S) decays recorded with the BABAR detector at the SLAC B factory. Identification of 
the B + — ► K + vv final state, with two neutrinos, requires the reconstruction of the companion B in 
the event. The companion B is reconstructed in the decay mode B~~ — > D°£~V X, which provides 
both high efficiency and good purity. The particles not used in the reconstruction of the companion 
B are compared with the signature expected for B + — > K + vV decays. Two candidates are found 
in the data with an expected background of 2.2 events. Under the assumption that all candidates 
are signal events, an upper limit on the branching fraction for B + — > K + vV of 9.4 x 10 -5 at 90% 
confidence level is determined. 
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1 Introduction 



The investigation of flavor-changing neutral current (FCNC) decays is of fundamental interest. In 
the Standard Model (SM) these decays are forbidden at tree level, and occur only in loop diagrams. 
As a result, their rates are highly suppressed. The SM prediction for the FCNC decay b — ► svV is 
nearly free from strong interaction effects and has very small theoretical uncertainty. An observation 
of this decay at a level significantly above the SM prediction would provide unambiguous evidence 
for new physics. 

Within the SM the decay b — > svV proceeds through W box diagrams and Z penguin diagrams. 
The expected branching fraction, summed over all neutrino species, is Q 

Bib ^sw) = (4.ll?;{j) x 1(T 5 . (1) 

At present it does not appear to be feasible to search for the inclusive decay b — > svV ; however, 
the decay B + — > K + vV is tractable.^ The expected branching fraction for B + — > K + vV , summed 
over all neutrino species, is |2j 

B(B + -» K + vV ) = (0.38ljj:Jl) x 10- 5 . (2) 

The best previous experimental limit is B{B+ -> K + uV ) < 2.4 x 10 -4 at 90% confidence level @. 

2 The dataset 

The data used in this analysis were collected with the BABAR detector, which is described else- 
where pi], at the PEP-II storage ring. The integrated luminosity used in this analysis is 50.7 fb -1 
recorded at the T(4S) resonance, corresponding to 56.3 x 10 6 BB events, and 6.4 fb -1 taken at en- 
ergies just below BB threshold. Simulated data samples for the processes e + e~ — > BB, e + e~~ — > qq 
(q = u, d, s or c) and e + e~ — > r + r _ , in quantities comparable to the data, are used to study 
backgrounds. A sample of 280 000 simulated B + B~ events with B + — > K + vv and the other S 
decaying generically have also been analyzed. The simulation of B + — > decay is based on 

the form factor model in Ref. 0. 

3 Analysis method 

The presence of two neutrinos in the final state makes the search for B + — > K^vv difficult, 
since no kinematic constraints can be applied to the signal B. The strategy adopted in this 
analysis is to reconstruct exclusively the decay of one of the B mesons in the event, referred 
to as the "tag" B, and to compare the remaining particle(s) in the event with the signature 
expected for the decay B + — > K + vV . The low multiplicity of the signal decay greatly reduces 
the combinatorial background in the tag reconstruction, allowing the use of decay modes that 
would not be sufficiently clean in other circumstances. These considerations lead to the use of 
the semileptonic decay B — > D°£~17 X for the reconstruction of the tag B. The X system is 
kinematically constrained to be either nothing or a low-momentum pion or photon from a higher 
mass charm state. The D° is reconstructed in the K~tt + , K~ir + Tr~ir + and K~it + it® modes. This 
method results in roughly 0.5% of B~ decays being reconstructed as tags. Note that particles from 

1 Charge conjugate modes are implied throughout. 
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the tag B that escape detection will not affect the sensitivity of the analysis to B + — > K + vV events; 
the reconstructed Di needs to be a correct, but not complete, subset of the particles produced in 
the tag B decay. The feed-down from higher-mass charm states often results in good tags in this 
sense, and thus in an enhanced tagging efficiency. 

The event selection proceeds as follows. Selected hadronic events are required to have an 
identified electron or muon with a momentum above 1.3 GeV/c in the T(4S) rest frame. The 
electron identification is based on quantities from the electromagnetic calorimeter (EMC), the 
ring-imaging Cherenkov detector (DIRC) and the gas (DCH) and silicon (SVT) tracking devices. 
The muon identification uses information from the instrumented flux return (IFR) in addition to 
the devices listed previously. Loose consistency requirements are placed on the charged particle 
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Figure 1: The candidate D invariant mass distributions are shown, from left to right, in the K~ir + , 
K~ir + ir~ir + and K~tt + tt° modes for data (points) and simulation (histogram), for events with no 
more than three charged tracks and less than 1 GeV of neutral energy not assigned to the tag B 
candidate. Events are required to have no more than three charged tracks not associated with the 
tag B in order to mimic the low multiplicity of the signal while maintaining adequate statistics in 
the plots. The off-resonance distributions have been scaled to the on-peak data luminosity. 



vertices for the D° and D°£~ candidates. The following kinematic requirements are imposed: 
p* D o > 0.5GeV/c, m D oe- > 3GeV/c 2 and —2.5 < cos 6 b,di < 1.1, where p* D0 is the momentum of 
the D° in the T(4S) frame, m^oi- is the mass of the D°£~ combination and 

. 2E B E m -m\-m 2 m 

cos B ,Dt = TP — i ■ (3) 

2 \PB I \PDt I 

Here Eb and \pb\ axe respectively the energy of and magnitude of the momentum of the B meson 
in the T(4S) frame. Eb is one half of the center-of-mass energy of the e + e~ initial state, and \pb\ 



is y (E B — m B ). The upper limit on cos 8b,D£ is 11 to account for resolution on the measurement 
(the signal cannot exceed 1). The lower limit is relaxed to increase efficiency for the feed-down 
from decays of the type B~ — ► D*°i~V and B~ — > D**°£~V. The requirement on cos 6b, di is the 
most important for restricting the kinematics of the D°£~ to be consistent with coming from a 
semileptonic B decay. In cases where more than one D°£~ candidate is reconstructed, the one with 
the smallest value of jcosf^^ | is used. The reconstructed D° invariant mass distributions are 
shown in Fig. |l[ 
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Once the tag B is selected, additional requirements are placed on the remaining particles in 
the event. There must be exactly one charged track in the event that is not part of the tag B, its 
charge must be opposite to that of the tag lepton, and it must satisfy the particle identification 
criteria for charged kaons, which are based on information from the DIRC and tracking system. 
The momentum spectrum, in the T(45) rest frame, for the kaon from B + — > K + uV decays peaks 
near the upper kinematic limit while the spectrum for background peaks at low momentum; the 
signal kaon candidate is thus required to satisfy p* K > 1.5GeV/c (see Fig. |2|). The angle 9* K g 
between the charged lepton and the signal kaon is isotropically distributed in signal events, since 
these particles originate from different B mesons, while the background from e + e~ — > qq and 
e + e~ — > t + t~ peaks forward and backward in this angle; we require —0.9 < cos9* Kg < 0.8. In 
addition to the above requirements on charged tracks, we use information from the EMC and 
IFR to limit additional neutral particles in the event. The B + — > K + uV signal leaves very little 
neutral energy in the detector and does not contain any neutral hadrons. We therefore require that 
the number of IFR clusters consistent with neutral hadrons (-/Vtfr) be zero, and that the energy 
deposited in the EMC, once the daughters from the D£ have been removed, (referred to as E\ e f t or 
remaining neutral energy) satisfies E\ e f t < 0.5 GeV (see Fig. |2|). 




|3500 - 




Remaining neutral energy (GeV) 



Figure 2: The distributions of p* K and E\ e ft for simulated signal and background events. Events 
with no more than three charged tracks and less than 1 GeV of neutral energy not assigned to the 
tag B candidate are used for the plot on the left whereas the neutral energy requirement is relaxed 
to be less than 3 GeV for the plot on the right. The generic MC distribution has been scaled to 
the on-peak data luminosity with an arbitrary scale factor applied to the signal MC distribution. 

The yields in the signal and sideband regions at each stage in the application of the selection 
criteria are given in Table [j] for the on-peak data and background Monte Carlo, along with the 
efficiency for the signal Monte Carlo. The distribution of events in the search plane defined by 
the variables^ E\ e f t and (uid — m %)/^ is shown in Fig. ||. The signal box is defined by the 

2 The quantities m% and a% are the mean and sigma from Gaussian fits to the D° invariant mass spectrum. 
Separate values are calculated for each D° decay mode in data and simulation. 
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Figure 3: The distribution of events in the (rnn — ™%)/p^d versus E\ e f t plane for on-peak data, 
generic BB and continuum Monte Carlo and signal Monte Carlo. In the generic Monte Carlo 
plot the circles show the contribution from BB events, the squares show the contribution from cc 
and the triangles show the contribution from uu/dd/s's. The MC has not been scaled to the data 
luminosity. 
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Table 1: The number of events passing the selection criteria for on-peak data, on-peak Monte Carlo 
contributions, off-peak data, off-peak Monte Carlo contributions and B + — > K + vV signal Monte 
Carlo efficiency. The number of events in the Monte Carlo sample are scaled to the equivalent 
luminosity in data. The values include the correction factors for tag efficiency, E\ e ^ and iVipR 
referred to in the text. 





On-peak (50.7 ib -1 ) 


Off-peak 


[6.4 ftr 1 ) 


signal MC 


Requirement 


data yield 


MC yield 


data yield 


MC yield 


effic -10 4 


Tag, no extra tracks 


8998 


8525.7 


415 


389.9 


34.3 


Kaon identification 


717 


707.4 


49 


46.8 


24.3 


COS0^ 


485 


486.2 


32 


25.0 


20.9 


* 

Pk 


101 


89.4 


7 


5.1 


14.2 


N WK 


79 


72.5 


6 


4.4 


12.0 


Eieft sideband 


34 


27.4 


3 


1.4 


0.2 


L> u mass sideband 


4 


7.1 


1 


0.8 


2.0 


Signal box 


2 


2.2 





0.3 


10.3 



requirements E\ e f t < 0.5 GeV and \m D — rrv^\ < 30""). The expected background from the Monte 
Carlo is 2.2 events. 

In order to minimize experimental bias, the signal region was hidden until the selection criteria 
were finalized. In order to evaluate how well the simulation describes the data, we define auxiliary 
samples. Two sideband regions are studied: the D° mass sideband, defined by the conditions 
\mr> — > 3cr^ and E\ e { t < 0.5 GeV, and a sideband where the additional neutral energy is 
required to be in the range 1.0 < E\ c { t < 2.5 GeV. The D° mass sideband contains incorrectly 
reconstructed B decays and continuum events whereas the E\ e f t sideband is sensitive to correctly 
reconstructed B tags where the other B leaves only a single detected charged track and substantial 
missing energy, often in the form of neutral hadrons. The event yields in these regions are also 
listed in Table [|. 

In addition to the sideband samples, we use "double-tagged" events, in which both B + and B~ 
mesons are reconstructed as B — > DiuX, to quantify the uncertainty in the efficiency of several 
of our signal criteria. We reconstruct double-tagged events by finding a suitable D°£~ candidate 
where the D° decays to K~n + , and then looking for a second D°£ + candidate in any of the accepted 
D° modes. No particle is assigned to more than one of the D£ candidates. In addition it is required 
that the event contain no charged tracks that are not assigned to a Di candidate. 

The reconstructed invariant masses of the D° and D° candidates in double-tagged events that 
satisfy the above criteria are shown in Fig. Q for data and for Monte Carlo. After subtracting 
combinatorial background, the number of double-tagged events satisfying the requirement that 
\rriD — < 3<jp for each D candidate is 148 ± 15 in data and 175 ± 16 in the Monte Carlo 
sample.^ The number of double-tagged events per fb _1 in the data is 0.85 ±0.11 times the rate 
in the simulation. This factor is roughly the square of the data/Monte Carlo efficiency ratio for 
the tag efficiency (including the requirement that there be no additional charged tracks associated 
with the tag - see the first entry in Table p]) • The signal efficiency is therefore corrected by a factor 
0.92 ± 0.06, where the uncertainty is taken as a systematic error. 

3 The number of events in the Monte Carlo sample has been scaled to the on-peak data luminosity. 
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Figure 4: Mass of the D candidate decaying to all three modes considered versus the D° — > ivT _ 7r + 
mass for events in which both B mesons are reconstructed in the D£ vX decay mode and no 
additional charged particles are recorded. The plot on the left (right) shows the results from the 
simulation (on-peak data). 
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Figure 5: The distribution of E\ e f t for "double-tagged" events where both B mesons are recon- 
structed in the Dt X decay mode and no additional charged particles are recorded. The plots on 
the left (right) show the distribution from simulation (on-peak data). 
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Table 2: A summary of the systematic uncertainties on B{B + — > K + vv ). 5e/e is the relative 
uncertainty on the overall efficiency. 



Quantity 


Se/e[%] 


_B.B-counting 


1.1 


Tagging efficiency 


6.0 


K selection 


2.0 


cos6* Kl 




Eleit 


4.3 


NlFR 


3.6 


K momentum 


1.8 



The double-tagged events also allow a study of how well the variables A?ifr and -Ei e f t are 
simulated. Figure || shows the distribution of the E\ e { t variable in the double-tagged events; the 
D° mass sidebands have been used to subtract the combinatorial background. The mean values of 
.Eieft in the data and simulation are 0.91 ±0.08 GeV and 0.84 ±0.07 GeV, respectively. The fraction 
of double-tagged events satisfying the requirement -/Vtfr = is 0.87 ± 0.03 in data and 0.93 ± 0.02 
in simulation. These comparisons are used to adjust the simulated signal efficiencies and assign 
systematic errors. 

Systematic uncertainties on the efficiency of selection criteria based on the total number of events 
with T(4S) mesons, tagging efficiency, K selection and momentum, -EWt and -/Vtfr have all been 
studied. The total relative uncertainty on the selection efficiency is found to be 5e/e = 8.7% where 
the tagging efficiency and -EWt contribute the largest uncertainties. The systematic uncertainties 
are summarised in Table |2[ 

4 Physics results 

The signal region was unblinded to reveal two events, consistent with the 2.2 events predicted with 
the simulation. The number of B + — > K + vV candidates in the data is thus compatible with the 
background expectation. For the purpose of setting an upper limit, each candidate is assumed to 
be signal. The Poisson upper limit for 2 events is 5.3. This upper limit must be modified to account 
for the uncertainty in the efficiency. Using the prescription advocated in [|| increases the upper 
limit to 5.4 events, from which we find 

B{B + -> K + vV ) < 9.4 x 10~ 5 (preliminary) (4) 

at 90% confidence level. 

The background at present appears to be mostly combinatorial, based for example on the lack 
of any D° peak in the continuum in Fig. |l]. Further refinements may enable this background to be 
suppressed; the combinatorial component of the background can also be subtracted in the future. 

The approach used in this analysis can easily be extended to B° — > K a vV and B — > K*uV as 
well as to B + — > t + u . 
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